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Abstract 

Using a simple model for nuclear GPDs, we study the role of the neutron contribution to nuclear 

"oo ' DVCS observables. As an example, we use the beam-spin asymmetry Af^^ measured in coherent 

O ' 

■ and incoherent DVCS on a wide range of nuclear targets in the HERMES and JLab kinematics. We 

(N : , 

^ , find that at small values of the momentum transfer t, Af^jj is dominated by the coherent-enriched 

contribution, which enhances Af^jj compared to the free proton asymmetry j4^^, Af^u{(f)) / A^jj{(j)) = 



■ 1.8 — 2.2. At large values of t, the nuclear asymmetry is dominated by the incoherent contribution 

and Af^jj/{(j))A^^jj{(j)) = 0.66 - 0.74. The deviation of Af^jj{(j))/Al^{(j)) from unity at large t is a 
result of the neutron contribution, which gives a possibility to constain neutron GPDs in incoherent 

■ 

^ ' nuclear DVCS. A similar trend is expected for other DVCS asymmetries. 
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I. INTRODUCTION 
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Hard exclusive reactions such as Deeply Virtual Compton Scattering (DVCS), 7*T 7T', 
and hard exclusive meson production (HEMP), 7*T MT', have emerged as indispensable 
s to access the microscopic (parton) properties of hadrons 

1^ . [13I . [l^ . In the above reactions, T and T' stand for any hadronic target (nucleon, 
pion, atomic nucleus); M denotes any meson. Note that the above reactions may also 
include transitions between different hadronic states such as e.g. N ^ A, p n, N 



Nn 



m 
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171] and production of pairs of mesons 



In the Bjorken limit (large Q^), the QCD factorization theorem for DVCS and HEMP 



on any hadronic target 



13 



141 ] states that corresponding scattering amplitudes factorize 
in convolution of perturbative (hard) coefficient functions with nonperturbative (soft) ma- 
trix elements, which are parameterized in terms of generalized parton distributions (GPDs). 
GPDs are universal (process-independent) functions that contain information on parton dis- 
tributions and correlations in hadrons and in matrix elements describing transitions between 
different hadrons (see above). 

In this paper, we consider DVCS on nuclear targets, 'y*A 7A, which gives an access 
to nuclear GPDs. We would like to single out the following three important roles of nuclear 
DVCS: 

• It gives information on nucleon GPDs, which is complimentary to that obtained in 
DVCS on the free proton; 

• It allows to study novel nuclear effects, which decouple from DIS and elastic scattering 
on nuclei; 

• It imposes stringent constraints on theoretical models attempting to give a covariant 
description of the nuclear structure. 

In this paper, we deal with the first point. In particular, we examine the role of the 
neutron contribution to nuclear DVCS asymmetries on a wide range of nuclei. This allows 
one to constrain neutron GPDs, which are not directly accessible. 

Nuclear DVCS opens possibilities to study novel nuclear effects, which seem to be pre- 
dominantly encoded in the real part of the DVCS scattering amplitude. It was speculated 
in the framework of the nuclear liquid drop model that the so-called nuclear Z)-term, which 



contributes to the real part of the nuclear DVCS amplitude, has a fast, non-trivial depen- 
dence on the atomic number A {A^^'^ vs. naively expected A^) lOL This observation was 
confirmed by an analysis of nuclear GPDs using the Walecka model 20| . In that analysis, the 
fast y4-dependence of nuclear GPDs comes from nuclear meson degrees of freedom. Hence, 
the measurement of DVCS observables sensitive to the real part of the DVCS amplitude 
gives a possibility to study non-nucleon (mesonic) degrees of freedom in nuclei. 

In the small Bjorken xb limit, a model for nuclear GPDs, which combines the model for 
nucleon GPDs based on the aligned-jet model with phenomenological parameterizations of 



usual nuclear PDFs, was suggested in 



21 



22| |. It was found that the ratio of the real parts 



of the nuclear to nucleon DVCS amplitudes has a very unexpected behavior as a function of 
Xb, which is very different from the corresponding ratio of the imaginary parts. The latter 
was found to be similar to the ratio of the nuclear to nucleon structure functions measured 
in inclusive DIS. This, again, hints that novel nuclear effects might be lurking in the real 
part of the nuclear DVCS amplitude. 

The third role of nuclear DVCS is related to the fact that nuclear GPDs, similarly to 
nucleon GPDs, should obey the fundamental properties of polynomiality and positivity. In 
order to achieve these properties, theoretical models used to build nuclear GPDs must give 
a covariant description of the nuclear structure, which imposes severe constraints on the 



nuclear models. This problem was discussed in relation to modeling deuteron GPDs in |23 |. 

The literature on nuclear DVCS and nuclear GPDs in not numerous and can be readily 
comprehensively overviewed. 

Originally, the formalism of deuteron GPDs was developed in 2^. The formalism of 
nuclear GPDs of any spin-0, spin-1/2 and spin-1 nuclei was presented in {2^. Assuming 
that nuclei are collections of free protons and neutrons, predictions for DVCS observables 
(asymmetries) were made. In particular, in accord with the earlier result of |26|, it was 
predicted that the nuclear DVCS beam-spin asymmetry is enhanced compared to the free 
proton asymmetry, A^jji^cj)) / A^^^cj)) ~ 5/3, for spin-0 and spin-1/2 nuclei. 

Up to now, the main theoretical approach to dynamical models of nuclear GPDs is the 
convolution approximation, which assumes that nuclear GPDs are given by convolution of 
unmodified or modified nucleon GPDs with the distribution of nucleons in the nuclear target. 
The latter distribution is obtained from the non-relativistic nuclear wave function. Within 



the convolution approximation, there were considered GPDs of such nuclei as deuterium 23|, 



3 



27 



208 



Pb 



28|^He |29|, l30(, ^He |3l|, [3^, ^^Ne and ^^Kr |26|, a wide range of nuclei from to 



20| (in that analysis, besides nucleons, meson degrees of freedom were also used in 



the convolution). 

While the convolution approximation is reliable for > 0.1, it is not applicable for small 
xb, where such coherent nuclear effects as nuclear shadowing and antishadowing become 
important. A model of nuclear GPDs for heav y n uclei, which takes into account nuclear 
shadowing and antishadowing, was proposed in [2l|, |22| (see also the discussion above). 

Another important aspect of nuclear DVCS, at least from the practical point of view, is 
the interplay between the coherent (the nucleus remains intact) and incoherent (the nucleus 
excites or breaks up) contributions to nuclear DVCS. This was studied in 26| and a general 
expression for nuclear DVCS asymmetries, which interpolates between the coherent and in- 
coherent regimes, was derived. It was predicted that for the coherent contribution, in the 
kinematics of the HERMES experiment, the ratio of the nuclear (^°Ne and ''^Kr) to the free 
proton beam-spin asymmetries is enhanced, Af^(0)/A^^(0) ^ 1.8. For the incoherent con- 
tribution, it was predicted that 74f^(0)/y4^^(0) = 1, provided that the neutron contribution 
to the nuclear DVCS amplitude is neglected. 

It is the main goal of the present work to go beyond this approximation and to study 
the role of the neutron contribution in coherent and incoherent nuclear DVCS observables 
(asymmetries). 



On the experimental side, initial 
HERMES collaboration at DESY 



measurements of nuclear DVCS were reported by the 



331] and more data on nuclear DVCS at HERMES is 



expected 3J]. The CLAS collaboration at Jefferson Lab recently reported a measurement 
of DVCS on deuterium with the aim to study the neutron GPDs 35|. It is planned that 
nuclear GPDs will be studied at Jefferson Lab at the present 6 GeV and the future 12 GeV 
energy of the electron beam. At high energies, nuclear GPDs will be studied at the LHC 
in ultraperipheral nucleus- nucleus collisions, see e.g. [361], and at the future Electron-Ion 
Collider. 

This paper is organized as follows. In Sect. [Tll we explain our model of nuclear GPDs. 
The interpolating formula between the coherent and incoherent regimes of nuclear DVCS 
is derived in Sect. IIIII Predictions for the nuclear beam-spin DVCS asymmetry in HER- 
MES and JLab kinematics, with an emphasis on the neutron contribution, are presented in 
Sect. HVl We summarize and discuss our results in Sect. |Vl 
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II. MODEL FOR NUCLEAR AND NUCLEON GPDS 



We use a simple model for nuclear GPDs that captures main features of the dependence 
of nuclear GPDs on the atomic number A and on the momentum transfer t. We assume that 
the nucleus consists of A uncorrelated nucleons: Z protons and N = A — Z neutrons 
see Fig. [H 



25|, 




FIG. 1: Schematic representation of nuclear quark GPDs. 



For simplicity, we shall consider spin-0 nuclei. In this case, there is only one leading-twist 
quark nuclear GPD, H'^, which can be expressed in terms of the free proton and neutron 
quark GPDs H'^ and E"^ as follows. 



FA{t) , (1) 



where -Fa(^) is the nuclear form factor normalized to unity; m^r is the nucleon mass; other 
variables are introduced below. Note that the GPDs H'^ and E"^ enter Eq. ([T]) in the com- 
bination that leads to the proper nuclear charge form factor 26 1. 

The Bjorken variable xa is defined with respect to the nuclear target. In the laboratory 
frame, we have 



Xa 



■Xb 



2uMa 2uAmN A 
where u is the photon energy; Ma is the mass of the nucleus. From the relations 



(2) 



it follows that 



-^ = A^^. (4) 

Next we find the relation between x and xjy. In the symmetric notation [81], the outgoing 
interacting quark carries the plus-momentum = {x + ^^)-P4 , see the left-hand side of 
Fig. [H On the other hand, can also be written as (see the right-hand side of Fig. [1]) 

k+ = {xn + e7v)P^ = {xn + ^n) (^jPl + ^) 

= ixN + ^N)(^jil + U)-U^Pl. (5) 

In this derivation, we used the assumption that = P^/^- Therefore, with help of Eq. (jlj), 
we find that 

In the forward limit, Eq. ([l]) reduces to the model for nuclear quark parton distribution 
functions (PDFs), 

qA{xA,Q^) = A[Z qpixB,Q^) + N qn{xB,Q^)] • (7) 
These nuclear PDFs satisfy the baryon number (total charge) and momentum sum rules, 

j dxA^eqqAixA,Q^) = j dxB^eq[Z qp{xB,Q^) + N qn{xB,Q^)] = Z , 



j dxA^XAqAixA.Q"^) = j dxB^XB 



Z N 



q ^ q 

Taking the first x-moment of the nuclear GPD weighted with quark charges, one obtains 
the nuclear electric form factor, 

F^it) = f dxY, eq Kix, U, Q\ t) = [ZF|(i) + NF^it)] Fa^I) , (9) 

where F^^{t) = Ff'"'{t) + t/(4m^)Ff '"(t) are the electric form factors of the proton and 
neutron expressed in terms of the corresponding Dirac and Pauli form factors. 

The fact that the right-hand side of Eq. ([9]) does not depend on C,a means that the first 
x-moment of H'^ satisfies polynomiality. An examination shows that higher x-moments of 
Hj^ do not satisfy polynomiality, even if the proton and neutron GPDs do. As we mentioned 
in the Introduction, it is an outstanding theoretical challenge to build a model of nuclear 
GPDs with the property of polynomiality. 



DVCS observables are expressed in terms of the so-called Compton form factors (CFFs), 
which are defined as nuclear GPDs convoluted with the corresponding hard scattering coef- 
ficients. For spin-0 nuclei, to the leading order in a^, the only OFF reads 



1 



1 



X - ^A + iO X + C 



iO 



6v 



dx 



N 



Z (n^/^ixN, Q^ t) + -J-^E'^/P{x^, Q\ t)^ 



+ N(^H'^/^{xN,^N,Q\t) 



4m% 



E^/^{xN,iN.Q\t) 



X FA{t) 



6v 



1 



XN-^N + i^ XN+^N-i^ 



t 



FA{t) . 



(10) 



An important corollary of Eq. ( |TOl) is that I-La scales as A^. 

In our analysis, for the nucleon CFFs Ti^'" and f^^'", we used results of the dual parame- 
terization of nucleon GPDs [37|, which gives a good description of the data on DVCS cross 



section and DVCS asymmetries on the proton target. In the modeling of the nucleon GPD 
E, we took Ju = Jd = 0. 



For the nuclear form factor FA(t), for ^He, we used the result of 38!]. For other nuclei, 
we used the parameterization of nuclear charge density distributions |33] (see Appendix for 
details) . 



III. COHERENT AND INCOHERENT NUCLEAR DVCS 

In the situation, when the recoiled nucleus is not detected, measurements of DVCS observ- 



ables with nuclear targets necessarily involve the coherent and incoherent contributions |26 |. 
The former contribution corresponds to the case when the nuclear target stays intact, and 
it dominates at small values of the momentum transfer t. The latter contribution corre- 
sponds to the case when the initial nucleus A transforms into the system of A — 1 spectator 
nucleons (bound or free) and one interacting nucleon, and it dominates at large t. In the 
approximation of closure over the final nuclear states, the exact structure of the final system 




FIG. 2: The coherent DVCS (left) and Bethe-Heitler (right) scattering amphtudes on a nucleus A. 
Only one of two possible BH amplitudes is shown. 

of A — 1 nucleons is not important. The coherent DVCS and BH amplitudes (one of two 
possible attachments of the real photon to the lepton lines is shown) are presented in Fig. [2l 
the incoherent DVCS and BH amplitudes are shown in Fig. [31 



FIG. 3: The incoherent DVCS and Bethe-Heitler scattering amplitudes. The initial nucleus A 
transforms into a final state containing A — 1 spectator nucleons (free or bound) and an interacting 
nucleon. 

In order to correctly sum the coherent and incoherent contributions to the eA — >■ e'yA 

□ 

cross section, let us schematically write the corresponding amplitude 




A 



A 



A-1 



A 




(11) 
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where A* represents the final state consisting of A nucleons (coherently scattered nucleus 
or any product of the nuclear dissociation); Jj represents the operator corresponding to the 
interaction with the nucleon i (one-particle operator); the summation runs over all nucleons 
of the target; A is the momentum transfer. Assuming that the states \A*) form a complete 
set, the cross section summed over the nuclear final states can be expressed in the following 
form, 



daA 
~dt 



AAA 

A* j i i,j 

A A 

= {A\ ^ Jj J. e''^-(--^-^)| A) + {A\ Jl MA) 

^ A{A - 1) {A\JI Jn e*^-(^^-^^ V) + A {N\jI J^|iV) 

where daN /dt is the scattering cross section on the bound nucleon; da^/dt corresponds to 
the quasi-free nucleon; t' = A/{A — l)t [40]. For the sake of the argument, we did not 
distinguish between protons and neutrons. Adopting the HERMES terminology, we shall 

. The second term is the 



call the first term in the last line of Eq. (fT2!) coherent- enriched 
incoherent contribution. 

The dependence of the coherent-enriched contribution on t is steep and is governed by 
the nuclear form factor squared F^{t'). Therefore, this contribution dominates the nuclear 
cross section at small t. The t-dependence of the incoherent contribution is much slower 
and is determined by the t-dependence of the cross section on quasi- free nucleons da^/dt. 
While this contribution is present at all t, it dominates the nuclear cross section at large t. 

Besides the t-dependence, the coherent-enriched and incoherent contributions have differ- 
ent ^-dependences. The coherent-enriched contribution scales as A{A — 1); the incoherent 
contribution scales as A. 

Let us now consider the case, when the recoiled nucleus is intact. In this case, 1^4*) = \A) 
in Eq. f|TT]) . and the expression for the eA — * e'-fA cross section becomes 

^ = A'Flit)^. (13) 
dt ^ dt ^ ^ 

In Eq. ( IT3|) . daA/dt is the genuine coherent nuclear scattering cross section, which scales as 
A^ and whose t-dependence is steep and is given by the nuclear form factor squared F\{t). 



Using Eq. (fT2!) . the full-fledged differential cross section for the eA e^yA reaction [iC 
can be written as a sum of the coherent-enriched and incoherent contributions, 



daA 



a^XAU 



A-l 



dxAdydtd<i) ^TtQ'^^lTe^ \ A 



^AixA^y) 



+ 



xb 
xa 



2 A 

E 

i=l 



%{xB,y) 



(14) 



where T4 is the amplitude for the coherent eA e'yA scattering; % are the amplitudes for 
quasi- free incoherent eA — > 67^4 scattering; the {A — I) /A factor originates from Eq. f|T2l) : 
the {xb/xaY factor is required for the incoherent contribution not to depend on A; is the 
angle between the lepton scattering and the production planes. 

It is important to note that the prefactor {A—l)/A corresponds to the DVCS amplitude 
squared and to the interference between the DVCS and Bethe-Heitler (BH) amplitudes. For 
the BH amplitude squared, {A — I) /A should be replaced by {Z — l)/Z. 

In Eq. ( IT4|) . in the laboratory frame, 

^XaMa ^XbTTIn 



y 



V 

'e 



(15) 



Q Q ' 

where E is the energy (momentum) of the incoming lepton. Note that the variables y, e and 
t are the same for nuclear and nucleon targets. 

For the comparison with the free nucleon case and with experiments, it is convenient to 
express a a as a function of xb, 

A 



daA 



a xbV 



A 



dxBdydtd(f) SnQ^VTT^ \ 



%i{xA,y) 



i=l 



%{xB,y) 



(16) 



For the BH amplitude squared, {A — 1)/A^ should be replaced by {Z — 1)/{ZA'^). 

For illustration, let us consider the DVCS contribution to Eq. fITBl) . In this case, |73ip oc 
which scales as [A'^FAit')]'^, see Eq. f ITU]) . Therefore, the first term in Eq. ([TBI) behaves 
as A'^F]^{t'). The second term has the t-dependence determined by the nucleon GPDs and 
scales as A. 

In the situation, when the recoiled nucleus is detected, the eA —>■ ejA cross section is 
purely coherent, and it reads 



daA 



a XBy 



dxBdydtdcf) 8ttQ^VTT~^ A'^ 



T^A{xA,y) 



(17) 
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IV. NUCLEAR DVCS ASYMMETRIES 



In this section, as an example of DVCS asymmetries, we consider the beam-spin nuclear 
DVCS asymmetry in the presence of the coherent and incoherent contributions, with an 
emphasis on the neutron contribution. We make predictions relevant for the HERMES and 
JLab kinematics. 



A. Coherent and incoherent contributions to DVCS asymmetries 

Expressions for nuclear DVCS asymmetries can be readily obtained from Eqs. f|T6|) and 
(fT7|) . In this work, we consider the beam-spin asymmetry, A^u, which is measured with the 
longitudinally-polarized lepton beam and the unpolarized target. 

The nuclear and nucleon amplitudes squared in Eqs. (fT6ll and (ITTll receive contributions 
from the DVCS and Bethe-Heitler (BH) scattering amplitudes and their interference, 

|T|'= |Tdvcs|'+|Tbh|'+X, (18) 

where I = '^vcs'^h + '^*h'^vcs- 



-'BH + -L + L/d""" 



The expression for the nuclear DVCS beam-spin asymmetry reads lOj 

AJ 

'-BU\' -rT+ I iDVCS 

where AJ = 1/2(J^=^ - J-^=-^) with A the helicity of the incoming lepton; all other contri 
butions correspond to the unpolarized beam. 

In the situation corresponding to Eq. ([TBI) , each term in Eq. (HM contains the coherent 
enriched and incoherent contributions, 

IT |2 _ ^ ~ ^ \q-A |2 I 7 \q-p |2 I Tvr \q-n |2 

I^bhI - ry ^2 I^bhI +^I^bhI +^^I^bhI ■> 



ZA^ 

^ A — 1 ^ 2 2 2 

I'^vcsl = ^3 I'^vcsl +^l'^vcsl +^l'^vcsl ■ (20) 



Expressions for the free nucleon contributions X^'"^, \T^l^\'^ and I'^vcsl^ terms of cos and 
sin (/)-harmonics are derived in [lOil. As a model of the nucleon GPDs, we used the results of 



the dual parameterization of nucleon GPDs with J„ = = 37 1. 



11 



Expressions for T^, I'^bhP I^^dvcsP spin-0 zero nuclei are the same as for the 
pion [41], after the replacement of the pion charge form factor by the nuclear one evaluated 
at t' = A/{A - l)t. 

In the case of the purely coherent scattering corresponding to Eq. f[T71) . the terms in 
Eq. (I2U]) should be replaced by the following expressions, 

X = —I^ 

i/T- |2 |2 

I^bhI - ^I^bhI , 

I'^DVCSP = -^\%DVCs\'^- (21) 

In the purely coherent case, the nuclear form factor is evaluated at the momentum transfer 
t. 

Using Eqs. (fT9l) . (!20l) and fl2T]) . one can qualitatively estimate the behavior of as 
a function of A and Z. Provided the |7BHp-term dominates the unpolarized cross section, 
the coherent-enriched contribution to Af^jj{(j)) scales as {A — 1)/{Z — 1). The purely coherent 
Af^jj{(j)) scales as A/Z. 

All expressions used in Eqs. fl20|) and fl2Tl) are collected in Appendix. 



B. Nuclear DVCS beam-spin asymmetry Am in HERMES kinematics 

In the measurement of nuclear DVCS at HERMES, the recoiled nucleus is not detected. 



3J]. This corresponds to the situ- 



but reconstructed using the missing mass technique 
ation, when one sums over all final nuclear states. This means that the nuclear beam-spin 
DVCS asymmetry, Af;^, receives contribution from the coherent-enriched and incoherent 
terms and is given by Eqs. (|T9l) and (1201) . 

We quantify our numerical predictions for Af^^ by considering the ratio of the nuclear to 
the free proton asymmetries, Af^^{(f)) / A^j^jj{(j)) . This ratio is presented in Fig. |4]as a function 



of t at an average HERMES kinematic point xb = 0.065 and = 1.7 GeV^ 3J|. The 
asymmetries are evaluated at = 90°. Different curves correspond to different nuclei: ^He, 
i^N, 20Ne, 84Kr and ^^iXe. 

At small values of t, when the nuclear asymmetries (cross sections) are dominated by the 
coherent-enriched contribution, Af^jj/{(j))A^j^jj{(f)) = 1.8 — 2.2, which is consistent with the 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

-t [GeV^] 

FIG. 4: The ratio of the nuclear to free proton beam-spin DVCS asymmetries, ^^f;((/>)/A^jy (</)), as 
a function of the momentum transfer t for He, N, Ne, Kr and Xe nuclei. The calculation is done 
atxB = 0.065, Q2 = 1.7 Q^y^ ^ and (/> = 90°. 



previous analyses 25|, |26|. The enhancement of y4fj;(0)/yl^^(0) above unity is the combi- 



natoric effect: Since the interference between the Bethe-Heitler and the DVCS amplitudes 
scales as Z{A — 1) and the Bethe-Heitler amplitude squared scales as Z{Z — 1), Af^^{(j)) 
scales as {A — 1) / {Z — 1). 

At large values of t, when the nuclear form factor eliminates the coherent-enriched term, 
^Lui4>) is given by the incoherent contribution, and Af^(0)/y4^^(0) < 1. 

The fact that y4ff;(0)/y4^^(0) < 1 is a result of the neutron contribution to Af^^{(f)), see 
Eq. (!20!) . First (this is effect is largest), the neutron contribution decreases the numerator of 
Af^jj^cj)), since Fin < 0, while Fip > 0. Second, the positive neutron contributions |7^jjp+T" 
(somewhat suppressed by the neutron electromagnetic form factors compared to the proton 
contribution) and I^dVcsP (similar to the proton contribution) increase the denominator of 
^Lui't')- The decrease of the numerator of Af^^{(f)) and the increase of the denominator work 
together to reduce Af;u{(j)) / A^^jj{(f)) significantly below unity at large t. 

Note that our present finding that A'^^i^cl)) / A^'j^jj^cf)) < 1 at large t does not contradict the 
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FIG. 5: The ratio of the nuclear to free proton beam-spin DVCS asymmetries, Af^^{(j))/A^^{(|)), as 
a function of the momentum transfer t for He, N, Ne, Kr and Xe nuclei. The calculation is done at 
xb = 0.065, = 1.7 GeV^ [3j] and = 90°. The left panel corresponds to the coherent-enriched 



contribution; the right panel - to the incoherent contribution. 



original analysis 



261] . In that work, is was predicted that Af^^{(j))/Al 



1 as t becomes 



large, if the neutron contribution to the nuclear asymmetry is neglected. In the present 
work, we went beyond this approximation and found that the neutron contribution is not 
negligible and leads to Af^^{(f)) / A^J^^{(f)) < 1. Therefore, studies of the incoherent contribution 
to nuclear DVCS asymmetries is a sensitive tool to constrain neutron GPDs. The CLAS 



35|. 



collaboration at Jefferson Lab explored this possibility using the deuterium target 

Note also that the neutron GPDs enter the model of nuclear GPDs, see Eq. ([T]). Hence, 
nuclear DVCS observables in the coherent regime also provide certain constraints for the neu- 
tron GPDs, albeit those constraints are less stringent and more model-dependent compared 
to the incoherent regime. 

By studying the t-dependence of the nuclear DVCS cross section, the HERMES analysis 
separated the coherent-enriched and incoherent contributions to Af^u{(f)). Our predictions 
for these two contributions are presented separately in Fig. [5l The left panel corresponds to 
the coherent-enriched contribution to Af^jj^cp), which was calculated keeping only first terms 
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FIG. 6: The ratio of the nuclear to free proton beam-spin DVCS asymmetries, Af;^{(j))/Al^ {(/)), 
as a function of A. The calculation is done at xb = 0.065, = 1.7 GeV^ and cp = 90°. 



in Eq. fl20l) . The right panel corresponds to the incoherent contribution calculated using last 
two terms in Eq. (120|) . 

In the left panel of Fig. \5[ the curve for ^He lies above the curves for other nuclei because 
the coherent-enriched contribution scales {A — 1)/{Z — 1). 

In the right panel of Fig. [5l the ratio Af^jj{(f)) / A^^jj{(f)) at small t is close to unity because the 
neutron contribution is suppressed by the small value of the neutron Dirac form factor 
As I Fin (if:) I increases with increasing the ratio y4^^(0)/A^(y(0) begins to progressively 
deviate from unity. 

Taking different t-slices of Fig. HI we can study the A-dependence of Af^^{(f)). Figure [6] 
presents Af^jj{(j)) / A^^jj{(j)) as a function of A at t = —0.018 GeV^ (upper set of points) and 
t = —0.2 GeV^ (lower set of points). These two values of t correspond to the average 



HERMES values 



34|. 



The interpretation of Fig. [6] is the same as for Fig. HI At small values of t, the coherent- 
enriched contribution dominates and Af^^{(j)) / A^J^^{(f)) > 1 due to the fact that Af^^{(j)) scales 
roughly as {A — 1)/{Z — 1). At large t, where only the incoherent contribution matters, 
^Lui'i') / ^^Lui'^) < ^ neutron contribution (see the discussion above). 

Results presented in Fig. [6] should be compared to the results of the HERMES analy- 
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SIS 



34| . At t = —0.018 GeV^, the agreement between our calculations (the upper set of 



points) and the HERMES data is excellent. For the nuclei of ^He, ^^N, ^ONe and ^"^Kr, 
Aj;^{(l))/Al^{(p) ^ 1.65. For the nucleus of ^^iXe, = 1.23, which is smaller 

than for other lighter nuclei because of the reduction of the coherent-enriched contribution 
by the nuclear form factor. 

At t = -0.2 GeV^, we predict that = 0.66-0.74, depending of the target 

nucleus. The experimental uncertainties of the HERMES data are too large and, in general, 
do not exclude the deviation of Af^jj{(f)) / A^^^{(j)) from unity, as we predict. 



C. Nuclear DVCS beam-spin asymmetry Am in Jefferson Lab kinematics 

There exists an exciting possibility to study purely coherent nuclear DVCS at Jefferson 
Lab using the BoNuS recoil detector. In particular, an experiment to study Am in coherent 
and incoherent DVCS on ^He has been proposed. Main advantages of the proposed experi- 
ment compared to HERMES is exclusivity of the measurement, which will allow to measure 
the purely coherent DVCS, and small projected errors due to high statistics, which will en- 
able one to unambiguously determine the magnitude of Am in the coherent and incoherent 
regimes. In addition, the proposed experiment might shed some light on the question of 
modifications of nucleon GPDs in nuclear medium. 

Figure [7] presents our predictions for the ratio of the coherent ^He to free proton beam- 
spin DVCS asymmetries, as a function of the momentum transfer t. The 
calculation corresponds to a typical point in the current JLab kinematics: E = 6 GeV, 
xb = 0.15 and = 1.5 GeV^. The ratio of the asymmetries is evaluated at = 90°. 

The behavior of A'^^^cf)) / A^j^jj^cl)) presented in Fig. [7] is similar to that in the left-hand side 
panel of Fig. O Since the purely coherent scales as A/Z, while the coherent-enriched 

contribution to scales as {A — 1)/[Z — 1), the curve in Fig. [7] lies lower than the 

corresponding curve in the left-hand side of Fig. [51 

In relation to incoherent DVCS on ^He, the proposed experiment at Jefferson Lab will 
measure the e^He epX reaction, i.e. DVCS on a quasi-free proton. In this case, the 
neutron contribution is absent and the ratio y4f^(0)/A^^(0) = 1, provided the bound proton 
in ^He is not modified. Therefore, this measurement will probe modifications of proton GPDs 
in ^He. 
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FIG. 7: The ratio of the coherent ^He to free proton beam-spin DVCS asymmetries, 
Af^jj{(f>)/A^jj{cl)), as a function of the momentum transfer t. The calculation corresponds to JLab 
kinematics, E = 6 GeV, xb = 0.15, = 1.5 GeV^, and was performed at (j) = 90°. 



V. CONCLUSIONS AND DISCUSSION 

Using a simple model for nuclear GPDs, we studied the role of the neutron contribution to 
nuclear DVCS observables. As an example, we used the beam-spin asymmetry Af^^ measured 
in coherent and incoherent DVCS on a wide range of nuclear targets. In our analysis, we 
considered the ^He, ^^N, ^''Ne, ^'^Kr and ^^^Xe nuclei in the HERMES kinematics and the 
^He nucleus in the JLab kinematics. 

We found that at small values of the momentum transfer t, Aj^^ is dominated by the 
coherent-enriched contribution, which scales as {A — 1)/{Z — 1). This enhances the nuclear 
A^^ compared to the free proton A^j;, Af^jj{(j)) / A^jj{(j)) = 1.8 — 2.2, in accord with the earlier 



predictions 



25 



26|. 



On the other hand, at large values of t, when the nuclear asymmetry is domi- 
nated by the incoherent contribution, Af^u{(j)) / A^^jj{(j)) is significantly smaller than unity: 
Af^u{(j))/A%{(j)) = 0.66 - 0.74, depending on the target nucleus. This deviation of 
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Af^^{(f)) / A^J^^{(f)) from unity is a result of the neutron contribution: The negative neutron 
contribution < 0) decreases the numerator of Af^jj and, at the same time, the positive 
neutron contribution I'T^hP + 2^" + I^^VcsP increases the denominator of Af^jj. Since the 
effect of the deviation of Af^^{(f)) / A^j^jj{(f)) from unity is so sizable, incoherent DVCS on nuclei 
gives a possibility to constain neutron GPDs. 

In this work, we considered one kind of DVCS observables, namely, the beam-spin asym- 
metry. We expect that for other DVCS asymmetries, such as e.g. for the beam-charge DVCS 
asymmetry, the ratio of the nuclear to the free proton asymmetries will be qualitatively sim- 
ilar to see |26|. 

All results presented in this work, data grids for the dual parameterization of the nu- 
cleon GPDs and FORTRAN codes for various DVCS asymmetries measured in DVCS 
on nucleons and nuclei can be found and downloaded from the author's webpage 
http : // www . j 1 ab . org/ ~ vguzey. 
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APPENDIX: INPUT FOR CALCULATION OF DVCS ASYMMETRIES 

In this Appendix, we collect all expressions used in our numerical analysis of the nuclear 
and proton DVCS beam-spin asymmetries, see Eqs. (fT9|) . (!20!) and (1211) . 

The interference, Bethe-Heitler and DVCS terms, which enter Eqs. ( IT9l) . ( !20l) and ( !2TI) . 
read Q 



IT |2 _ _^ „DVCS /A IN 
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where we have kept only twist-two terms and neglected the gluon GPDs. In Eq. ( 1A.1|) . Vi{4>) 
and 7^2(0) are lepton propagators; the plus-sign in front of the interference term corresponds 
to electrons, while the minus-sign is for positrons; is the angle between the lepton and 
production planes; the coefficients Cq ^j^p, s{ ^^p, '^o^np called harmonics. 

When Eq. (lA.ip is applied to the coherent-enriched contribution, it should be evaluated 
with X = xa and the corresponding nuclear harmonics (see below). When Eq. (lA.ip is used 
to calculate the incoherent contribution, it should be evaluated with with x = xb and with 
the free proton and neutron harmonics (see below). 



1. Nuclear part 



Expressions for the cos(f) and sin 0-harmonics of a spin-0 zero nucleus are the same as 
for the pion case 4l|] after the replacement of the pion GPD and the electric form factor 
by their nuclear counterparts (one has also divide the pion harmonics involving GPDs by 
the factor of x due to a different normalization of the interference and DVCS terms used 



m 



4l|). The required harmonics read 
„ t 



^0,unp 

''l,unp 
-.2: 

'l,unp 

.BH 
-'0,unp 



Q 
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ZFA{t) ^cHa 



-8K{2-2y + y^) ZFA{t) ^cHa 
%K\y{2-y) ZFA{t) '^niHA, 

{{2-yy + y\l + er) 



t 

+ 4(1 - Xa) + {AxA + e^) — 
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+ 2 [4(1 - y)(3 + 2e') + y\2 - e')] - Ax\{2 - yf{2 + e')— )Z'Fl{t) 
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^2,unp 

DVCS 
OjUnp 



-^K{2-y) (^2xA + e'-4xl^ 
32K^x^, ^ 



t 



2{2-2y + y') \nA\\ 



(A.2) 



where K is the so-called kinematic i^-factor lo|; A is the incoming lepton helicity. 

The nuclear form factor Fa entering Eq. ( ]A.2[) is evaluated at t' = A/{A — l)t for the 
coherent-enriched contribution and at t for the purely coherent case. 
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For ^He, the nuclear form factor is parameterized as 



FA{t) = (1 - (ahf) e 



(A.3) 



where a = 0.316 fm and b = 0.681 fm 



For other nuclei used in this paper, the nuclear form factor is defined as 

Jo v t 



(A.4) 



where pAi'f') is the nuclear charge density distribution taken in the following form 391]. 
Nitrogen {A = U, Z = 7): 



w = -0.18, z = 0.505, 



2.57, po = 0.0127908 



pA{r) = po 



1 + e(^-'=)/^ ' 



Neon {A = 20, Z = 10): 



Po 
PA{r) 



Krypton (A = 8A, Z = 36): 



z 

Po 
PA{r) 



Xenon {A = 131, Z = 54): 



0.571, c = 2.805, 



0.00767524 , 



Po 



1 _|_ (,{r-c)/z 



0.496, c = 4.83, 
0.00191897, 



Po 



1 + e('"-'=)/^ ' 



(A.5) 



(A.6) 



(A.7) 



w 
c 



PAir) 



0.3749 , z = 2.6776 , 
5.3376, Po = 0.00112617, 

Po: 



1 + 



I _|_ g(r2-c2)/22 



(A.8) 
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2. Proton part 



ffressions for the required cos0 and sin 0-harmonics for the proton target are derived 
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where 
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Equations (1A.9P and (lA.lOp involve proton Compton form factors (CFFs) Tip and £p and 



electromagnetic form factors Fip and F2p. For the CFFs, we used the dual parameterization 
with Ju = Jd = ^ [s^l- The proton electromagnetic form factors are parameterized in the 
following form [l^ 

k 

F2p{t) = YZ ~^D{t), 

Gn{t) = j^, (A.ll) 
where kp is the proton anomalous magnetic moment, kp = 1.79; my = 0.84 GeV. More 



elaborate parameterizations of the nucleon elastic form factors are possible, see e.g. 4^, but 
Eq. flA.lip is sufficiently accurate for our purposes. 

3. Neutron part 

Expressions for the cos0 and sin 0-harmonics for the neutron case are readily obtained 
from Eqs. (]A.9I) and (lA.lOl) by replacing the proton CFFs and electromagnetic form factors 
by their neutron counterparts. The neutron CFFs are obtained from the proton ones by 
exchanging e„ ^ in the DVCS amplitude. 

The neutron electromagnetic form factors are parameterized in the following form [1( 

t k 

Fln{t) = —Goit), 

F2n{t) = -^^Goit), (A.12) 
where kn is the neutron anomalous magnetic moment, /c„ = —1.91. 
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